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ABSTRACT 

We run adiabatic A'^-body/hydrodynamical simulations of isolated self-gravitating gas 
clouds to test whether conformal gravity, an alternative theory to General Relativity, 
is able to explain the properties of X-ray galaxy clusters without resorting to dark 
matter. We show that the gas clouds rapidly reach equilibrium with a density profile 
which is well fit by a /3-model whose normalization and slope are in approximate 
agreement with observations. However, conformal gravity fails to yield the observed 
thermal properties of the gas cloud: (i) the mean temperature is at least an order of 
magnitude larger than observed; (ii) the temperature profiles increase with the square 
of the distance from the cluster center, in clear disagreement with real X-ray clusters. 
These results depend on a gravitational potential whose parameters reproduce the 
velocity rotation curves of spiral galaxies. However, this parametrization stands on an 
arbitrarily chosen conformal factor. It remains to be seen whether a different conformal 
factor, specified by a spontaneous breaking of the conformal symmetry, can reconcile 
this theory with observations. 

Key vifords: gravitation - methods: A^-body simulations - galaxies: clusters: general 
- dark matter - X-ray: galaxies: clusters. 



1 INTRODUCTION 

On the scale of individual galaxies and larger scales. Gen- 
eral Relativity requires large amounts of dark matter to 
describe the dynamics of cosmic structure. Moreover, the 
late-time acceleration of the Hubble expansion implies the 
existence of a cosmological constant, a special case of a dark 
energy fluid which is suggest ed by more sophisticated mod- 
els fsee ICopeland et al.ll2006l . for a review). In principle, we 
can avoid the dark matter and dark energy solutions to the 
puzzles posed by the astrophysical data by adopting an al- 
ternative theory of gravity, which reduces to General Rel- 
ativity in the appropriate limit. Independently of the dark 
matter and dark energy problems, a modification of General 
Relativity is also highly desirable if we ultimately wish to 
unify gravity with the other fundamental interactions. 

Possible modifications of the Einstein-Hilbert action 
proposed in the literature are, among others, (i) the in- 
troduction of additional scalar and/or vector fields (e.g., 
iFuiii fc Maeda 2003: Bekenstein 2006): (ii) the assumption 
of arbitrary functions f(R) of the Ricci scalar R (e.| 
ICapozziello fc Francaviglial l2008l : iNoiiri fc Odintsovl 12008 
(iii) the introduction of additional dimensions to the four di- 
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mensions of the General Relativity spacetime manifold (e.g., 
iMaartendliooi ). 

A different approach was s uggested by iMannh eimI 

l|l990l ) who revived Weyl's theory (wevllllQia [1919, 192^ 
as a possible candidate to solve the dark matter and dark en- 
ergy problems. When the geometry is kept Riemannian, with 
a null covariant derivative of the metric tensor, we can obtain 
a milder version of Weyl's gravity, known as conformal grav- 
ity. In this theory, we impose a local conformal invariance on 
the gravitational field action in the curved four-dimensional 
spacetime. The Einstein-Hilbert Lagrangian density for the 
gravitational field is chosen on the requirement that the the- 
ory of gravity is a second-order derivative theory. In con- 
formal gravity, the Lagrangian density is chosen on the 
principle of local conformal symmetry which is uniquely 
satisfied by the action Iw = ~ct J d'^Xyf—gC^vK.xC^"'^^ , 
where Cfj,uK\ is the Weyl tensor, a is a coupling constant, 
and g is the determinant of the metric tensor g^i^. Con- 
formal symmetry is garanteed by the invariance of the 
Weyl tensor to local conformal transformations guvix) —* 
Q^{x)gfj,,y{x), where Q^{x) is an arbitrary conformal fac- 
tor that can be speci fied by a spontaneous breaking of th e 
conformal invariance lEderv et al.ll2006l : lMannheimll2008bh . 
The theory of gravity implied by the action Iw is a fourth- 
order derivative theory. The vacuum exterior solution for 
a static and spherically symmetric spacetime contains the 
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Schwarzschild solution (|Mannheim fc Kazanas|[l989l ). The 
weak-field lim it is consistent w ith the Solar System obser- 
vational data (lMannheimll2007|). unlike cla i med in previous 
investigations (Barabash fc Shtanovl Il999l : iFlanaganI l2006l : 
iBarabash fc P vatkovska"2007''l . 

Mannheim ( 1993, , 1997 ) shows that conformal gravity 
can reproduce the rotation curves of disc galaxies without 
dark matter. Moreover, conformal gravity can explain the 
current accelerated expansion of the universe without re- 
sorting to a fine-tun ed cosmological constant or to the exis - 
tence of dark energy ()Mannheimil200ll . l2003l : IVarieschj20oi ) . 
Unlike the standard theory, where the universe starts accel- 
erating at redshift z < 1, conformal gravity predicts that 
the universe is accelerating at all times. Therefore, obser- 
vatio nal data probing the Hubble plot at ve ry high redshift 
(e.g., iNavia et al.ll2008l : I Wei fc ZhariS |2008| ) can be a deci- 
sive test. 

More recently, conformal theory has been proposed as 
a v alid candidate for building a theory of quantum grav- 
ity (|Mannheiml l2008al ) : in fact, theories based on fourth- 
order derivative equations of motion have had the long- 
lasting problem of sufferin g from the presence of ghosts. 
iBender fc MannheimI (|200^) have recently shown that this 
erroneous belief is the result of considering the canonical 
conjugates p of the generic dynamical variables q, when the 
q's are real, to be Hermitian operators; however, this as- 
sumption is incorrect, and when the non-Hermiticity prop- 
erty of the Hamiltonian of these higher-order quantum field 
theories is correctly taken into account, the states with neg- 
ative norm disappear. 

From an astrophysical perspective, however, the form of 
conformal gravity that has been proposed in the literature 
currently has two main shortcomings: the abundance of light 
elements, and the gravitational lensing phenomenology. 

Conformal gravity nicely avoids the requirement of an 
initial Big Bang singularity, but it still predicts an early 
universe sufficiently dense and hot to ignite the light ele- 
ment nucleosynthesis. However, conformal gravity predicts 
a too slow initial expansion rate. This rate favours the de- 
struction of most of the deuterium produced in the early 
universe (|Knox fc Kosowskvl[l993l : lElizondo fc Yepeslll994 ') 
and poses conformal gravity in serious difficulties compared 
to the standard Big Bang nucleosynthesis. Conformal grav- 
ity necessarily requires astrophysical processes for the pro- 
duction of the deuterium currently observed, for example 
neutron radi ative capture on proton s in the atmospheres of 
active stars (IMullan fc Linskvl[l999l) or gamma-ray bursts 
l|lnoue et al.ll2003l '). However, the processes investigated to 
date do not seem to be efficient enough. For example, sig- 
nificant pr oduction of deuterium in t he Galaxy seems to be 
ruled out (jProdanovic fc FieldsllioO^ ) . 

The second open issue is the conformal gravity pre- 
diction of gravitational lensing. Early inv estigations of 
gravitational lensing in conformal gravity (jWalkeil 1 19941 : 
[Ederv fc Paraniapc 1998, 1999; Ederv 19^ show that, in 
the weak-field limit, the deflection angle due to a point mass 
M is Aa = 4GM/c^r — -yr, where r is the radius of the pho- 
ton closest approach; Aa contains the additional term "/r 
when compared to the General Relativity result. The con- 
stant 7 has to be positive to fit the galaxy rotation curves, 
thus implying a re pulsive effect in g ravitational lensing. It 
was later realized (|Ederv et al.ll200lh that the geodesies of 



photons are independent of the conformal factor (x) and 
one can choose an appropriate conformal factor and a radial 
coordinate transformation to yield attractive geodesies for 
both massive and massless particles. However, in the strong- 
field limit, the l ight defl ection might still be divergent or 
even impossible (|Pireauxi,2004a. .b,) . 

Until these open questions are completely settled, con- 
formal gravity cannot yet be considered ruled out by obser- 
vations. Moreover, conformal invariance plays a crucial role 
in elementary particle physics and a viable theory of grav- 
ity that includes this property can at least be suggestive of 
a relevant route towards the uniflcation of the fundamental 
interactions. 

From the astrophysical point of view, conformal grav- 
ity can be further tested by investigating the formation of 
cosmic structure. To date, nobody has yet explored how the 
large-scale structure forms in conformal gravity. If struc- 
tures form by gravitational instability, as in the standard 
theory, the development of a cosmological perturbation the- 
ory, which is still lacking, becomes inevitable. This theory 
would enable the comparison of conformal gravity with the 
spectrum of the Cosmic Microwave Background anisotropies 
and would provide initial conditions for the simulation of the 
structure evolution into the non-linear regime. 

Before building up such a theory, however, it is useful 
to check whether conformal gravity is able to reproduce the 
equilibrium configuration of cos mic structure s, other than 
galaxies, without dark matter. iHornd (|2006l ) has already 
shown that, if we interpret the observational data of the in- 
tracluster medium of X-ray clusters by assuming hydrostatic 
equilibrium, conformal gravity requires a factor of ten less 
baryonic mass than inferred from the X-ray surface bright- 
ness measures. Here, we extend this analysis by performing 
hydrodynamical simulations of self-gravitating gas clouds. 
We modify a standard Tree-|-SPH code to perform numeri- 
cal simulations of self-gravitating systems in conformal grav- 
ity. The numerical tool we create is extremely relevant if we 
eventually wish to investigate the formation of the large- 
scale structure, because we will massively need to resort to 
numerical integrations when the evolution of the density per- 
turbations reaches the non-linear regime. 

In Section [2] we review the basic steps leading to the 
gravitational potential of a static point source in confor- 
mal g ravity. In Section |31 inspired by the analysis of iHornj 
(120061 ) ■ we compute the gravitational potential energy of a 
spherical system, and in SectionUwe use this result to com- 
pute the expected mean temperature and the temperature 
profile of the intracluster medium. In Section (5] we derive 
the same results with hydrodynamical simulations of self- 
gravitating gas clouds. 



2 CONFORMAL GRAVITY 

Conformal gravity is a theory of gravity based on the action 
Iw = ~ct J d*x^—gCf_i^KxC^'^'^^ . For any finite, positive, 
non- vanishing, continuous, real function Q,^{x) of the four 
spacetime coordinates x, this theory is invariant for any con- 
formal transformation g^i^ix) Q,'^{x)g^v{x), because the 
Weyl tensor is invariant for these transformations. The vari- 
ational principle applied to the action leads to fourth-order 
field equations, whose solution for the vacuum exterior to a 



© 0000 RAS, MNRAS 000, 000-000 



X-ray clusters of galaxies in conformal gravity 3 



static point source leads to the line element 



ds 



2<^(r) 



1 + 



dr — r {d9 + sin 



where c is the speed of light, 



(r) 
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3/37) 



3^ , 7 k 2 
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and /3, 7 and k are three integration constants. The arbi- 
trary function Q?[x) can be specified by a mechanism which 
breaks the conformal symmetry. iMannheim fc Kazanad 
l|l989l) implicitly assume Q?(x) — 1. 

According to this assumption, when 7 = fe = 0, the 
metric Q reduces to the usual Schwarzschild metric with 
P — GM/c? , M the gravitational mass of the point source 
and G the gravitational constant. The term corre- 
sponds to a cosmological solution which is conformal to a 
Robertson- Walker background; therefore, k can be chosen 
small enough that it can be neglected on the scales of galax- 
ies and galaxy clusters, and we will not consider this term 
hereafter. Moreover, the rotation velocities of spiral galaxies 
suggest a value of 7 sufRciently small that the terms propor- 
tional to p-y oc 7/c^ can be safely ignored ()Mannheimyl99S) . 
Equation ((2)1 thus reduces to 



<t>{r) 



1 + 1, 



(3) 



[Mannhei^ |T993) uses the observed rotation curves of eleven 
spiral galaxies of widely different luminosities to constrain 



7. It tur ns out that 7 also 
mass M. iMannheimI l(l99i1 ) 



suggests the parametrization 
7 = 70 + Mjt, where 70 = 3.06 x 10~^* m""' and 7, = 
5.42 X 10"^^ m~^ ^0^; 70 and 7* should be two universal 
physical constants. In the following, we adop t a more conv e- 
nient form of the potential ^ suggested bv lHorn3 (|2006l ): 



GM GM GMo 



(4) 



where Mq = (7o/7,)M0 = 5.6 x lO^^M© and _Ro = 
{2GMQ/'y,c^y^'^ = 24 kpc. 



3 GRAVITATIONAL POTENTIAL ENERGY 
OF A SPHERICAL SYSTEM 

The gravitational potential of a point source (equation |4)| 

(j}ps{r) = 0]v(r-) + (j)CM{r) + (t>cc{r) , (5) 

with <j}N{r) = —GM/r the usual Newtonian potential, 
4>CM{r) = GMr/Ro, and (f>cc{r) = GMqt/R^, generates 
the acceleration 

, , ^ , GM GM GMo 
gps{r) = -V(j>j,s = -—2 ^2 (6) 

The last term is a constant acceleration independent of 
the mass of the source of the gravitational field. The phys- 
ical origin of this term is very subtle. Unlike in Newtonian 
gravity, the presence of the linear term cf)CM(r) clearly pre- 
vents us from neglecting the contribution of distant objects 
to the acceleration of any body in the universe. According to 



IMannheimI l| 19971 , 1 20061 ). in conformal gravity the net effect 
of all the mass in the universe is to contribute a constant 
acceler ation in addition to the acceleration due to a local 
source. IMannheimI shows that, with an appropriate 

coordinate transformation allowed by the conformal invari- 
ance, the mass in the universe exactly generates the addi- 
tional linear potential (jjccir), that we name the curvature 
potential. 

This argument implies that, to compute the gravita- 
tional potential energy of an extended source, we need 
to treat the curvature potential cf>cc{j') and the potential 
<i>N{r) + (f>cM{r), originated by the local field source, sepa- 
rately. For the time being, let us consider this latter gravi- 
tational potential alone. 

Consider an extended source. The gravitational poten- 
tial energy of two of its mass elements with separation r is 
dWAB = dmAdrriBfir), where /(r) is the part of the po- 
tential (j){r) that depends on r only. When we sum over all 
the mass elements, we get the total gravitational potential 
energy W = (1/2) J dniA J dmBf{r), namely 



(7) 



where p(r) is the mass density distribution of the source 
and <l){r) — 0jv(r) + ipcMijc) is the gravitational potential 
generated by its self-gravity alone. 

The total potential energy W can also be derived with a 
different argument which provides an alternative expression 
to equation ((T]). For a system of A'^ particles with position 
and mass rrii, each feeling a force F(ri), we can define the 
virial ^ , • F(ri); the force is F(ri) = —rui rij<i>{^ij) 
where Vij = Vi — Vj. Now, mi\7rij4>{rij) = — mj Vr^ . ^(rji) 
for Newton's third law, and we can write 



^r. ■F(rO = -Vm 



+ > rUi > Tj ■ Vr 



(8) 



since Vr^^ ^(rij) — when j — i. 

Homogeneous functions of order A of M variables 
(xi, • ■ ■ , xm) are defined by the relation 

/(axi, • ■ • , qxm) = Q^/(xi, • ■ • ,xm) (9) 

for any non-null a; they satisfy Euler's theorem 

Xj • Vi:;/(xi, • • • ,XAf) = A/(xi, . . . ,xa/) . (10) 



E 



In the potential (j) = (j)N + 4>cm, 0iv and (f)cM are homoge- 
neous functions of order Ajv = — 1 and Xcm ~ 1 respectively. 
By applying Euler's theorem for AI — 1, we can thus write 



E 



r^ ■ F(rO 



Wn - WcM 



mi[-~(j)Nirij) + 4>cM{rij)] 



(11) 

Now, W = Wn + WcM and we obtain W = 2Wn — J2i ' 
F{ri). In the continuous limit 



W = 2Wn + / p(r)r • V0(r)d^r 



(12) 
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Figure 1. Left panel: Total gravitational potential energy of a spherical system with a power-law density profile p cc r radius a = 1 
Mpc, and mass M = 2 X lO^^ M© (solid line), 4 X 10^^ Mq (dotted line), 6 X lO" Mq (dashed line), 8 X lO^^ M0 (dot-dashed Une), 
10^'^ Mq (dot-dot-dot-dashed line); see Section |3] for details. Right panel: The mean gas temperature of the systems shown in the left 
panel according to the virial theorem; see Section |4] for details. 



To compute the gravitational po tential 4)(r) o f a self- 
gravitating spherical system, we follow iHornel (|2006l l and we 
first consider a homogeneous spherical shell of density p, ra- 
dius R and mass m — AivpR^dR. At a generic point in space 
of coordinate r, each mass element Sm = pR^ sin ddRd9d(p 
of the shell generates the potential S<j)(r) given by equation 
((5|. By setting the coordinate system such that r — (0, 0, r) 
without loss of generality, the element Sm has coordinates 
i?(cos tfi sin 6, sin ip sin 9, cos 9) and generates the potential at 
r 



= GpR^ sin edRd9dip ( + 



(13) 



where = R^ + — 2rRcos9. By integrating over d9d<p 
we find the potential of the shell 

_ , „ / -* + 7i;(fi + «) 

*-"^"i-+^(*-) 

The gravitational potential of a self-gravitating sphere 
is thus 



(14) 



G 



loir) 



+ 



Rl 



3r 



+ r/o(r) 



-£;„i(r) + £i(r) 



(15) 



where lii{r) = A-n p{x)x"'^'^dx and 
An J^°° p{x)x"'^^dx, and the acceleration is 



E4r) = 



G 



Mr) 



+ 



Rl 



Hr) 
3r2 



-I- 



L(r)-7oCr) 



.(16) 



The total gravitational potential of a sphere can now 
be computed with either equation ([7]) or equation (|12|l . For a 
power-law density profile p(r) = po{r/a)~" of a system trun- 
cated at radius a with total mass M = 47rpoa;^/(3 — a), with 



Q / 3, we find 7o(r) = M{r/a)^-", h{r) = [(3 - q)/(5 
a)]Ma\r/af-", E-i{r) = [(3 - q)/(2 - a)](M/a)[l ■ 
(r/a)2-«], and £i(r) = Ma(3 - a)/(4 - q)[1 - (r/a)*"" 
and the gravitational potential energy is 



W = 



here 



Wn = 



Wn 
WnII 



(ro) 
hia)] ■ 



2(5 - 2a)(9 - la) 
3(7-2a)(5-Q) 



GM^ 3-Q 
a 5 — 2q: 



(17) 



(18) 



is the potential energy in Newtonian gravity, that is recov- 
ered in the limit Ro ^ 00. 

To compute the contribution Wcurv of the curvature ac- 
celeration — GMo/Rq to be added to M^, we can use either 
equation (|12|l or equation ([7]) without the factor (1/2), be- 
cause the curvature mass Mq is not accelerated by the mass 
p(r)d\. We find 



Wcurv = GMMq 



i?g4- 



= -Wn 



(ro) 



Mo 5 -2a 



RoJ M 4- 



(19) 



The left panel of Figure [T] shows the total gravitational 
potential energy 



WNia)[l~h{a)] + W,u,.ia) 



(20) 



of a spherical system as a function of a and its total mass 
M. When a 2.5,3.5,4 and 5, Wtot ±00. In Figure [1] 
the pole of Wn at a = 2.5 seems to have different properties 
than the other poles, but it is only a graphic artifact. In fact, 
the contribution of conformal gravity to Wtot is proportional 



to (a/Ro) 



10 , and it always dominates the newtonian 



Wjv: we must have \a — 2.5| <C 10^'' to see Wtot ±cx3. 

Finally, the ratio between the curvature potential and 
the conformal gravity potential due to the sphere self-gravity 
is ~ Mo/M; therefore, VFcurv is negligible when Ad 3> Mq, as 
in our case, where Mo/M ~ 10~^, unless, of course, a — > 4. 
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Figure 2. Upper left panel: Fluctuations of the total energy. Upper right panel: Evolution of the kinetic (dotted line), thermal (dashed 
line), gravitational potential (dot-dashed line), and total (solid line) energy. Lower left panel: Evolution of the mass- weighted mean X-ray 
temperature. Lower right panel: Evolution of the virial ratio. 



4 X-RAY CLUSTERS 

A system of A'^ particles of position experiencing the force 
Fi is in virial equilibrium when it satisfies the virial the- 
orem relation 2K -I- ^ ri • Fi = 0, where K is the total 
kinetic energy of the system. We saw above (equation llip 
that ^ Ti ■ Fi — 2Wn — Wtot in conformal gravity. Note 
that we replaced Wcm with Wtot to include the curvature 
potential. 

Observations indicate that in X-ray clusters the total 
mass of galaxies is roughly ten times s maller than the tota l 
mass of the intracluster medium (e.g., IWhite et al.lll993 ). 
We can therefore model X-ray clusters as individual clouds 
of hot gas with no galaxies and, of course, with no dark 
matter. This approximation is partly inappropriate for cD 
clusters, where the contribution of the cD galaxy to the grav- 
itational potential in the cluster center is relevant; however, 
cD clusters are ^ 20% of the population of nearby clusters 
l|Rood fc Sastrvlll97ll ) and our analysis should be valid for 
most clusters. 

If the gas bulk flows and turbulent motions are neg- 
ligible, the total kinetic energy is K = (3/2)(7 — 1)U = 
{3/2)Mk{T)/^imp, where U = [M/i-y - l)]k{T) / iinip is the 
internal energy of the gas of density p(r), M is the gas total 



mass, (T) = j p{r)T{r)(fir / J p{r)d'^r is its mass-weighted 
mean temperature, k the Boltzmann constant, mp the pro- 
ton mass, fi the mean molecular weight, and 7 the adiabatic 
index. The virial theorerrQ thus yields 

k{T) = -^{2W^-W,^,) (21) 

Figure [1] shows the mean temperatures k{T} of spherical 
systems with a power-law density profile as a function of a 
and its total mass AI. 

Consider an isothermal sphere with a = 2, radius a — 1 
Mpc, and mass Al = W^'^A4q. Newtonian gravity predicts 
k{T) = 0.09 keV. Conformal gravity predicts k{T) = 58.21 
keV, namely a temperature (a/Ro)^ ~ 10^ larger. To ex- 
plain a more typical temperature of, say, k{T) ~ 6 keV 
without resorting to dark matter, the gas mass should be 
Af ~ W^^AIq, a factor of ten lower than the mass of the in- 
tracluster medium measured from the X-ray surface bright- 
ness. 

This result agrees with the claim of lHornd (|2006l ) that 
the acceleration provided by the conformal gravity potential 

^ The virial theorem can also be generalized in f{R) gravities 
llBohmer et al.ll2008h . 
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is too intense. We now show that the disagreement with 
observations extends to the temperature profile, besides its 
normalization. 

To compute the temperature profile of a self-gravitating 
spherical system in equilibrium, we need to solve the Boltz- 
mann equation coupled to the Poisson equation. In confor- 
mal gravity, the Poisson equation V ^(/) = A-kGp is replaced 
by a fourth order equation V*</> cx p (|Mannheim fc Kazanad 
1 19941 ). Following this approach is not a trivial task, either 
analytically or numerically. In the next section, we resort to 
a smoothed-particle-hydrodynamical simulation to obtain a 
self-consistent solution to this problem. 

In this section, to illustrate qualitatively the expected 
result, we can consider the first moment of the Boltzmann 
equation, namely the Euler equation, which, in hydrostatic 
equilibrium, reduces to the relation Vp = — pV<^, where p 
is the gas pressure. Provided that the density profile p(r) is 
known, for an ideal gas with equation of state p = pkT / p.mp, 
this equation immediately integrates to 



fcr(r) 

piUp 



1 

p{r) 



p{r)\/(t){r)dr + Ao 



(22) 



where Aq is an integration constant. At large radii, real 
clusters usually have the intracluster medium density profile 
which decreases more rapidly than the temperature profile. 
This constraint requires that the term Ao/p{r) should dis- 
appear at large radii and implies Aq = 0. For a power-law 
density profile, we find 



T(r) = 



where 



TNir) 
4(3- 



1 + 



(2 -a) 



+ 



f-V 

\RoJ 
GMo 



(5- 



Ro 1 



fcTjv(r) _ GM 



/imp 



2a a 



IT (23) 



(24) 



Note that for an isothermal sphere (a = 2) in Newtonian 
gravity, this equation yields UTm / P-rrip = GM/2a. This tem- 
perature disagrees with the correct virial theorem result 
kT^/pmp = GM/3a (equation I21|) . This discrepancy is a 
simple consequence of the inconsistent approach we used to 
obtain equation (|23|) . However, we can use equation (|23|) to 
predict that, at large radii, the temperature profile increases 
at least as r . 

We thus conclude that the values of the conformal phys- 
ical constants Rq = 24 kpc and Mq = 5.6 x 10^° Mq imphed 
by the galaxy rotation curves provide X-ray clusters with 
average temperatures a factor of ten too large, and rising 
temperature profiles, clearly at odds with real clusters. 



5 NUMERICAL SOLUTION 

To solve properly the coupled Boltzmann and Poisson equa- 
tions and predict the correct X-ray properties of a galaxy 
cluster in conformal gravity, we run hydrodynamical simu- 
lations of isolated spheres of gas with a Smoothed-Particle- 
Hydrodynamic (SPH) code. We use a m odified version of the 
pubhcly availa ble GADGET-1.1 code (iSpringel et al.ll200ll : 
ISpringel|[2005l ). The original version of GADGET-1.1 is a 
Tree-|-SPH code which integrates the equations of motion 
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Figure 3. Evolution of the size of the gas cloud: From bottom to 
top, lines are for radii containing 10% to 90% of the total cloud 



of TV particles which interact gravitationally. The A'^ parti- 
cles are a discrete representation of either a coUisionless or 
a coUisional fluid; traditional simulations contain both col- 
lisionless and coUisional particles which represent the dark 
matter and the baryonic fluids, respectively. 

In our simulations, we only require the presence of 
the baryonic fluid and we only have coUisional particles. 
Moreover, in Newtonian gravity, the potential between 
two particles decreases with the interparticle distance, and 
GADGET-1.1 arranges the particles in a tree structure in 
order to average the contribution of distant particles to the 
local acceleration. This technique substantially decreases the 
time required to compute the acceleration of individual par- 
ticles. Unfortunately, in conformal gravity this technique can 
not be applied, because the gravitational potential contains 
a term that increases linearly with the interparticle distance, 
and the contribution of distant particles to the gravitational 
potential must be computed individually. We will therefore 
use GADGET-1.1 as a direct TV-body integrator. The com- 
putational cost thus increases with rather than TV In TV. 
In appendix [X] we describe in detail the modification we 
introduced in GADGET-1.1 to compute the interparticle 
forces in conformal gravity. 

We simulate the evolution of isolated spherical clouds of 
gas with vacuum boundary conditions. We neglect radiative 
cooling. At the end of this section we discuss how simulating 
the gas physics more realistically can actually be relevant to 
find a way to alleviate the discrepancy between conformal 
gravity and the observed thermal properties of X-ray clus- 
ters. 

A gas cloud is represented by TV = 4096 gas particles. 
The chosen TV is substantially smaller than state-of-the-art 
simulations of galaxy clusters and is due to the TV^ scaling of 
the computational cost. We do not actually need a larger TV, 
because the aim of our simulations is to test the general vi- 
ability of conformal gravity rather than the detailed physics 
of the intracluster medium, which would require higher spa- 
tial and mass resolutions. In fact, the current form of the 
conformal gravitational potential yields clusters so different 
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Figure 4. Mass density (left panel) and temperature (right panel) profiles of the final configuration of the gas cloud. The solid line in 
the left panel is the best-fit /3-model profile. The straight line in the right panel shows the T oc slope. 
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Figure 5. X-ray surface brightness (SB) (left panel) and emission weighted (EW) temperature (right panel) maps of the final configuration 
of the gas cloud. Each panel shows three orthogonal projections. Each map is 2 Mpc on a side. 



from real clusters that running costly simulations is unwar- 
ranted. 

The gas cloud initially has a /3-model density profile 



Pir) = po 



ro 



-313/2 



(25) 



and null thermal and kinetic energies. Each gas particle has 
gravitational softening length e = 0.1 kpc. Each gas parti- 
cle also has an individual SPH softening parameter h used 
to estimate the thermodynamical quantities, as mentioned 
below; the sphere with radius h centered on the particle 
contains a fixed number Nt of neighbor particles. We set 
A^f, — 40 in our simulations. 



In the following, we show the result of a typical sim- 
ulation: a cloud of total mass M = 2.7 x 10^^ M0, and 
initial values of the /3-model /3 = 0.636, ro — 134 kpc, and 
po = 1.52 X 10~^^ g cm~^. The upper left panel of Figure 
[5] shows that the total energy is conserved to better then 
0.1%. The evolution of the gravitational potential, kinetic, 
thermal and total energies (upper right panel) shows that 
in half a billion years the system reaches equilibrium and 
remains confined; in fact, 90% of the mass is within 500 kpc 
(Figure[3]). It follows that in conformal gravity a cloud of gas 
can easily remain confined without the need of dark matter. 

The bottom right panel of Figure [2] shows the evolution 
of the virial ratio: the sum K of the thermal and kinetic 
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energy over the total gravitational potential energy Wtot ■ At 
equilibrium we must have 2K + ^ • Fi = 0, where ^ • 
F» = 2Wn - VKtot (SectionBl). Therefore K/Wtot = (1/2) - 
(Wjv/VFtot)- From guation 1201 we see that Wtot is positive 
and Wtot ^ |Wjv| because [a/Ro)'^ ^ 1; therefore, unlike 
Newtonian gravity where K/Wtot ~ —1/2 at equilibrium, 
here the virial ratio K/Wtot is slightly larger than 1/2. 

In the Tree+SPH code each gas particles of mass rrii and 
density pi carries the electron number density riei, and the 
internal energy Ui . By assuming an optically thin gas of pri- 
mordial cosmic abundance X = 0.76 and Y — 0.24, we can 
estimate the particle temperature T;; we also estimate the 
ion number density nm + nuai that we use below. The bot- 
tom left panel of Figure [5]shows that the final mass- weighted 
mean temperature of the system is 90 keV, as expected from 
the argument provided in Section O 

Figure 13] shows the density and temperature profiles of 
the final configuration. As in real clusters, the /3— model is 
still an excellent approximation to the final density profile. 
The best fit yields fJ = 1.65, ro = 142 kpc and po = 1-6 x 
lO"'^^ g cm~^. These parameters roughly agree with typical 
values of real clusters, altho ugh /3 is ~ 2 — 3 times larger 
than typical values (see, e.g.. iMarkevitch et al.|[l999t ). 

The result in strong disagreement with observations is 
the rising temperature profile shown in the right panel of 
Figure 13] At large radii, we see that T oc r^, as anticipated 
in Section [4] 

For completeness, we show the surface brightness and 
the temperature maps in Figure [5] Each map is an equally 
spaced A^p x Np grid, with Np = 128, corresponding to a 
length resolution ~ 15.6 kpc. In the Tree-|-SPH code, each 
gas particle has a smoothing length hi and the thermo- 
dynamical quantities it carries are distributed within the 
sphere of radius hi according to the compact kernel W{r; hi), 
which has the same functional form of the gravitational ker- 
nel (see AppendixjXI , where r is the distance to the particle 
center. The X-ray surface brightness Sjk on the grid point 
{j, k} is 

^ ^ ^ ?^e^(nH^ -I- nHci)MTi)'WidVi (26) 

where dp is the pixel area, the sum runs over all the particles, 
and Wi cx J W{x)dl is the weight proportional to the fraction 
of the particle volume dVi = rrii/pi which contributes to 
the grid point {j, fc}. For each particle, the weights Wk are 
normalized to satisfy the relation ^ lUfc = 1 where the sum 
is now over the grid points within the particle circle. When 
hi is so small that the circle contains no grid point, the 
particle quantity is fully assigned to the closest grid point. 
We us e the cooling function A(r) of ISutherland fc Dopital 
l| 19931 ). The total X-ray luminosity is 

Lx=^ ueiinm + nH„)A(rOdVi . (27) 

Given the large equilibrium temperature, the cluster shown 
in Figure [5] has Lx ~ 3.86 x 10*'^ erg s~^ at equilibrium, 
an order of magnitude larger than the most luminous X-ray 
clusters. The temperature whose map is shown in Figured] 
is the emission-weighted temperature 

_ X) "ei('^Hi + nnBi)f^{Ti)TiWiAVi 



This map dramatically shows how the thermal properties of 
the intracluster medium totally disagree with real clusters. 

At this point, we cannot draw our conclusions with- 
out the following relevant cautionary tale. The results we 
show here derive from simulations that treat the gas adia- 
batically. In other words, we neglect radiative cooling and 
the gas heating processes due to astrophysical sources, as 
supernovae explosions, energy injection from active galac- 
tic nuclei, or galactic and stellar winds. These heating pro- 
cesses act at different times and with different effectiveness, 
depending on the detailed history of star and galaxy for- 
mation. This history is currently totally unexplored in con- 
formal gravity, and the detailed mechanisms of the heat- 
ing processes them selves are still poorly understood (e.g., 
iBorgani et al.l[2008l ). Therefore, the appropriate inclusion of 
these cooling and heating processes in our simulations is be- 
yond the illustrative purpose of this paper. 

However, an appropriate interplay between cooling and 
heating processes could in principle provide a possible route 
to reconcile conformal gravity with the properties of X-ray 
clusters. In fact, since our results show that, for the typi- 
cal gas density of real systems, the gas temperature is an 
order of magnitude larger than observed, radiative cooling 
would be very efficient in conformal clusters and the entire 
intracluster medium could cool in less than 1 Gyr. Effective 
heating processes would thus be necessary to reheat the intr- 
acluster medium to the observed X-ray temperatures. These 
heating processes must be more efficient in the cluster cen- 
ter than in the cluster outskirts, because the cooling time 
increases with the clustrocentric radius. However, although 
we cannot exclude that the various complex physical pro- 
cesses could eventually conspire to reconcile conformal grav- 
ity with observations, it is plausible that, in order to provide 
intracluster media with the observed thermal properties, the 
cooling and heating processes should be severely fine-tuned 
during the formation and evolution of clusters and cluster 
galaxies. 



6 CONCLUSION 

Conformal gravity can explain the rotation curves of disk 
galaxies and the current accelerated expansion of the uni- 
verse without resorting to dark matter and dark energy. We 
have modified a Tree-|-SPH code to run hydrodynamical sim- 
ulations of isolated X-ray galaxy clusters to show that con- 
formal gravity does not share the same success on the scales 
of clusters. 

These simulations confirm our simple analytic estimates 
that show that gas clouds with mass ~ 10^^ - 10^^ Mq, 
which are typical values of the total mass of the hot gas 
present in real clusters, remain confined with an equilibrium 
mean temperature ~ 10— 100 keV, ten times larger than the 
observed temperatures; more dramatically, because of the 
presence of a linear term in the gravitational potential, at 
large clustrocentric radius r, the gas temperature increases 
with r^, rather than decreasing as in real systems. 

Our analysis totally neglects radiative cooling and gas 
heating from astrophysical sources, as supernovae or active 
galactic nuclei. The interplay between these processes can in 
principle provide a way to reconcile conformal gravity with 
observations. It is however unclear if and how much these 
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processes should be fine-tuned to provide X-ray clusters in 
agreement with observations. 

In addition to this topic, we can see two more open 
issues whose solution might also reconcile conformal gravity 
with observations: 

(i) In conformal gravity, all the matter in the universe 
is expected to affect the local dynamics. The net effect is 
to contribute a constant inward acceleration —GMq/Ro in 
addition to the gravitational acceleration generated by lo- 
cal sources. Because we included this constant acceleration 
in our simulations, we could neglect the rest of the uni- 
verse and impose vacuum boundary conditions. It might 
be possible that assimilating the gravitational influence of 
the nearby matter surrounding the X-ray cluster in the con- 
stant "universe" acceleration ~GMq/Ro is inappropriate: 
in fact, nearby external matter might decrease the gravi- 
tational attraction of the interior matter and hopefully re- 
duce the thermal energy of the gas. To appropriately inves- 
tigate this effect, we should simulate the dynamics of large- 
scale structure within a full cosmological context. However, 
this task is not trivial just because the gravitational field is 
highly non-local. This investigation would also benefit from 
the implementation into the numerical simulation of the yet 
unavailable theory of structure formation. 

(ii) The gravitationl potential we implemented in our sim- 
ulations derives from a metric where the conformal invari- 
ance is broken by an arbitrary choice of the conformal factor 
Q^{x). It is unclear whether this choice provides a coordi- 
nate system whose physics describes the real world or it is 
an artifact of the reference frame. It also remains to be seen 
whether a spontaneous breaking of the conformal invariance, 
in th eories where matte r and gravity are conformally cou- 
pled l|Ederv et al.ll2006D . can provide a metric, and thus a 
gravitational potential, where the observed thermal proper- 
ties of the intracluster medium can be reproduced. 

In Section [T] we mentioned that the nucleosynthesis of 
light elements and the phenomenology of gravitational lens- 
ing are two open issues that need to be solved before ac- 
cepting conformal gravity as a viable alternative theory of 
gravity and cosmology. Here, we have shown that the ther- 
modynamics of X-ray clusters poses a third challenge to this 
theory. 



ACKNOWLEDGMENTS 

We thank Volker Springel for making public his superb 
Tree-|-SPH code and the referee whose careful report 
prompted us to clarify some aspects of our results. We ac- 
knowledge vivid discussions with Jorg Colberg, Keith Horne 
and Philip Mannheim on very early versions of this work. 
Support from the PRIN2006 grant "Costituenti fondamen- 
tali deirUniverso" of the Italian Ministry of University and 
Scientific Research and from the INFN grant PD51 is grate- 
fully acknowledged. 



Bekenstein J. D. 2006, Contemp. Phys., 47, 387 
Bender C. M., Mannheim P. D. 2008, PRL, 100, 110402 
Bohmer C. G., Harko T., Lobo F. S. N. 2008, JCAP, 03, 
024 

Borgani S., Diaferio A., Dolag K., Schindler S. 2008, Sp. 

Sci. Rev., 134, 269 
Capozziello S., Francaviglia M. 2008, Gen. Relativ. and 

Grav., 40, 357 

Copeland E. J., Sami M., Tsujikawa S. 2006, Int. J. Mod. 

Phys. D, 15, 1753 
Edery A. 1999, PRL, 83, 3990 

Edery A., Fabbri L., Paranjape M.B. 2006, Class. Quant. 
Grav., 23, 6409 

Edery A., Methot A.A., Paranjape M.B. 2001, Gen. Rela- 
tiv. and Grav., 33, 2075 

Edery A., Paranjape M.B. 1998, PRD, 58, 4011 

Edery A., Paranjape M.B. 1999, Gen. Relativ. and Grav., 
31, 1031 

Elizondo D., Yepes G. 1994, ApJ, 428, 17 
Flanagan E. E. 2006, PRD, 74, 023002 
Fujii Y., Maeda K. 2003, The scalar-tensor theory of grav- 
itation, Cambridge University Press 
Horne K. 2006, MNRAS, 369, 1667 

Inoue S., Iwamoto N., Orito M., Terasawa M. 2003, ApJ, 
595, 294 

Knox L., Kosowsky A. 1993, |arXiv:astro-ph/9311006 
Maartens R. 2004, Living Rev. ReL, 7, 7 



Mannheim 
Mannheim 
Mannheim 
Mannheim 
Mannheim 
Mannheim 

340 
Mannheim 
Mannheim 
Mannheim 
Mannheim 
Mannheim 

Grav., 26, 



D. 1990, Gen. Relativ. and Grav., 22, 289 

D. 1993, ApJ, 419, 150 

D. 1997, ApJ, 479, 659 

D. 2001, ApJ, 561, 1 

D. 2003, Int. J. Mod. Phys. D, 12, 893 

D. 2006, Prog, in Part, and Nucl. Phys., 56, 



D. 
D. 
D. 
D. 
D 
337 



2007, PRD, 75, 124006 

2008a, arXiv:0707.2283 

2008b, arXiv:0809.1200 

Kazanas D., 1989, ApJ, 342, 635 

. Kazanas D., 1994, Gen. Relativ. and 



Forman W. R., Sarazin C. L. 



H. 2008, 



REFERENCES 



Barabash O. V., PyatkovskaH. P. 2007, arXiv:0709.1044vl 
Barabash O. V., Shtanov Yu. V. 1999, PRD, 60, 064008 



Markevitch M., Vikhlinin A 

1999, ApJ, 527, 545 
MuUan D. J., Linsky J. L. 1999, ApJ, 511, 502 
Navia C. E., Augusto C. R. A., Tsui K. 
larXiv:0807 ":Q590l 

Nojiri S., Odintsov S. D. 2008, arXiv:0807.0685j 
Pireaux S. 2004a, Class. Quant. Grav., 21, 1897 
Pireaux S. 2004b, Class. Quant. Grav., 21, 4317 
Prodanovic T., Fields B. D. 2003, ApJ, 597, 48 
Rood H. J., Sastry G. N. 1971, PASP, 83, 313 
Springel V. 2005, MNRAS, 364, 1105 
Springel V., Yoshida N., White S. D. M. 2001, 

tronomy, 6, 51 
Sutherland R. S., Dopita M. A. 1993, ApJS, 88, 
Varieschi G. U. 2008, arXiv:0809.4729 
Walker M.A. 1994, ApJ, 430, 463 
Wei H., Zhang S. N. 2008, arXiv: 0808.2240l 
Weyl H. 1918, Math. Z., 2, 384 
Weyl H. 1919, Ann. Phys. Lpz., 59, 101 
Weyl H. 1920, Phys. Z., 21, 649 

White S. D. M., Navarro J. F., Evrard A. E., Frenk C. S. 
1993, Nature, 366, 429 



New As- 



253 



© 0000 RAS, MNRAS 000, 000-000 



10 Diaferio & Ostorero 



APPENDIX A: THE SOFTENING KERNEL 

The publicly availa ble code GAD GET- 1.1 (|Springel et alj 
I2OOII : ISprineelllioOsh is a Tree+SPH code which integrates 
the equations of motion of N particles which interact grav- 
itationally. The particles are a discrete representation of ei- 
ther a coUisionless or a colhsional fluid. Here, we consider 
only adiabatic processes and neglect the possibility of radia- 
tive cooling of the colhsional fluid. 

The only modification to the code we need is the com- 
putation of the gravitational potential (f> and its correspond- 
ing acceleration. At position r, A'^ particles of mass m; at 
position Xi generates the potential 



where yi = \r — Xi|. For point sources, gN(y) ~ 1/j/ and 
gcAiiy) = gcc{y) = y- The acceleration a(r) = —V(f){r) is 



where g]si{y)y = dgN/dy, and analogously for gcM and gcc- 
GADGET- 1.1 treats particles as extended spherical ob- 
jects with mass m and density profile p{r) — mW{r;h) 
where 



where x = r/h. W{r;h) is a spline kernel which avoids un- 
realistic divergences of the Newtonian acceleration for arbi- 
trary small particle separations. 

The spline kernel implies that each particle is not a 
point source of gravitational potential. Rather, its gravita- 
tional potential is correctly computed as an extended source 
of density W{r; h). According to equation (fT5)l we find, with 
u = y/h, 




N 





1 - 6a;2 + 62:^ O^x <\ 
2{l-xf \i^x<l 
x^l , 





gcM{y;h) = h < 




^ « < i 



i < U < 1 
U > 1 , 




gcc = y, and g^c = 
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